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Abstract 

Propagation, transmission and reflection properties of linearly po- 
larized plane waves and arbitrarily short electromagnetic pulses in one- 
dimensional dispersionless dielectric media possessing an arbitrary space- 
time dependence of the refractive index are studied by using a two-component, 
highly symmetric version of Maxwell's equations. The use of any slow 
varying amplitude approximation is avoided. Transfer matrices of sharp 
nonstationary interfaces are calculated explicitly, together with the am- 
plitudes of all secondary waves produced in the scattering. Time-varying 
multilayer structures and spatiotemporal lenses in various configurations 
are investigated analytically and numerically in a unified approach. Sev- 
eral new effects are reported, such as pulse compression, broadening and 
spectral manipulation of pulses by a spatiotemporal lens, and the closure 
of the forbidden frequency gaps with the subsequent opening of wavenum- 
ber bandgaps in a generalized Bragg reflector. 

1 Introduction 

Propagation of electromagnetic waves through stationary inhomogeneous media 
is a subject that has been extensively studied in the past 

The physics of moving dielectric media [3l [H O [H [7] is attracting renewed 
attention, mainly due to recent works by Visser [8] and Leonhardt and Piwnicki 
[9]. Particularly important is the formal equivalence between the equations of 
general relativity and the equations of light propagation in arbitrarily moving 
dielectric media, which has been fully recognized in the early literature fiO\. 

Much less attention, however, has been devoted to the case of light prop- 
agation through an inhomogeneous medium which possesses a time-dependent 
refractive index, and in particular light scattering by nonstationary interfaces 
between two media with different refractive indices. This must be distinguished 
from the above case of a moving medium, in that the nonstationary interfaces 
may possess arbitrary velocity in the range < < oo. This well-known 
fact [llj is not in conflict with special relativity because the medium itself is 
immobile, and only the interfaces of the transition regions of the refractive index 
change in time. This apparent motion of the interfaces is due to the time- varying 
nature of the refractive index, since it is not connected to the real motion of any 
physical particle or wave, and is not restricted by any limiting velocity, so that 
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the function n{z^t) that describes the spatiotemporal variations of the refractive 
index is completely arbitrary. In particular Lorentz transformations, which of 
course play a crucial role in the physics of moving media, are not applicable 
to the physics of media with a time- varying refractive index. For example, the 
field amplitudes of the waves generated inside the latter media are not found via 
Lorentz boosts, and the constitutive relations between the dielectric displace- 
ment and the electric field are fixed and not subject to Lorentz transformations, 
in contrast with the physics of moving media, described, for instance, in Refs. 

The first attempt to develop a theory for the very special case of plane 
wave propagation in a homogeneous medium with a sudden step-like temporal 
variation of refractive index has been given by Morgenthaler in 1958 [12j, but 
only after more than a decade this work has been noticed and other aspects of 
propagation in the same system have been explored by Felsen and Whitman in 
1970 [13j and by Fante in 1971 [14J. A more recent review of the salient features 
of the optics of nonstationary media, with particular emphasis on plasma physics 
applications, has been given by Shvartsburg in 2005 [15j. 

The purpose of the present paper is to study the effect of an arbitrary space- 
and time-dependence of refractive index n{z^t) on the propagation of electro- 
magnetic waves and pulses in media without dispersion. We shall provide a 
unified theoretical approach, with which we are able to study the scattering of 
light by nonstationary interfaces moving at arbitrary velocities. Light propaga- 
tion through media with a small number of nonstationary interfaces has been 
investigated by a variety of experimental techniques [19l [20l [211 EH [24] , as 
discussed further in the final discussion and conclusions section. Our results 
here include these cases but also address the behavior that could be achieved 
by further extension to an arbitrary number of well-controlled nonstationary 
interfaces. Our aim is to demonstrate the potential of nonstationary dielectrics 
for light control and manipulation. 

The paper is organized as follows. 

In Section [2] we construct the basic theoretical tools that are necessary for 
the analysis and the understanding of the building blocks of spatiotemporal 
dielectric structures. First of all, in subsection 12. II our master equation [Eq.([2])] 
is derived from Maxwell's equations, and its geometrical interpretation in terms 
of forward- and backward-moving fields is given. After that, in subsection 12.21 
the interface transfer matrix for the general case is calculated explicitly, while 
in subsection 12.31 the expression for free propagation when the fields are far 
from interfaces is derived, and a criterion for defining sharp interfaces is found. 
In subsection 12.41 we use special orthogonal variables to introduce plane wave 
expansion of arbitrary fields, which allows us to find the transfer propagation 
matrix in Fourier space. 

In Section [3] we describe three fundamental applications of the concept ex- 
pressed in Section[2l In subsection l3. li the refiection and transmission coefficients 
of several important configurations of nonstationary interfaces are derived, and 
numerical simulations of pulse scattering from interfaces are provided in order to 
confirm the theoretical expressions. In subsection 13.21 we discuss nonstationary 
multilayer structures, and provide the analogue of quarter wavelength condi- 
tion for a generalized nonstationary Bragg refiector, and the transformation of 
the bandstructure and the bandgaps under rotations of the (z, ct)-plane are ex- 
amined. In subsection 13.31 we deal with a special device, the spatiotemporal 
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analogue of a lens. We show, including numerical simulations, how such a lens 
can be used to perform pulse compression or broadening, with the consequent 
manipulation of pulse frequency and wave number. 

The results presented in Sections [2] and [3] are obtained by solving the full 
Maxwell equations. We derive in Section [H an exact equation to describe the 
evolution of the pulse envelope. We consider how this equation is modified if we 
introduce a slowly varying envelope approximation (SVEA) in time. Numerical 
simulations are then presented to show that the SVEA introduces errors in the 
detail of pulse propagation for the problems we are considering, due to the 
neglect of effects associated with sharp spatiotemporal interfaces. Finally we 
summarize our conclusions in Section O 

2 Analytical Results 
2.1 Main Equation 

Let us start by assuming that we have an electromagnetic wave, with its electric 
and magnetic fields E = [E{z^ t), 0, 0]^ and B = [0, B{z^ linearly polarized 

along the directions x and y respectively. The x-component of the electric 
field and the ^-component of the magnetic field depend on the variable z only, 
because in conditions of normal incidence we assume that any change of the time- 
dependent refractive index occurs along the propagation direction z. The linear 
polarization Pl of the medium will be written in the form Pl = x(^7^)E = 
[n(z, t)'^ — 1]E, where x is the linear susceptibility of the non-magnetic medium 
and n{z^t) is the linear space-time dependent refractive index, which is assumed 
for simplicity to be independent of the frequency co and to be real. This model 
for the susceptibility is obviously simplified, but for the purposes of the present 
paper the frequency dispersion of the refractive index is not essential and will 
be disregarded for sake of clarity. Here and in the following we choose to use the 
very convenient Heaviside-Lorentz units system [3 , in which E, B and Pl are 
measured in the same physical units (V/m), and there are no extra multiplying 
factors 47r, eo and jhq in Maxwell's equations. Maxwell's equations for the scalar 
quantities E{z^t) and B{z,t) can be written in the following form (c is the speed 
of light in vacuum) : 

dt{n^E)^cd^B = {), dtB ^ cd,E = 0, (1) 

where the two divergence Maxwell's equations are automatically satisfied due 
to the chosen polarization and the normal incidence assumption. 

It is possible to cast Eqs.([T]) in a different form, which turns out to be 
much more convenient for the analysis of time-varying photonic crystals and 
lenses, using the fact that the problem under consideration is (l+l)-dimensional. 
Defining local right- and left-moving fields, = {nE ± B)/2^ we can write 
Eqs.([T]) in the following compact form: 

* = 0, (2) 

where ^ = F~]^ (with the superscript T we indicate vector transposition), 
and 

= ( J ; ) , = ( ? J ) = ( J ) ' (3) 



d + (7x{dhin)cF] 
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are the identity matrix and the first and third (real) Pauh matrices. Moreover, 
= (ctq ±(Tj)/2 are the Pauh projectors {j = {x, z}), and d is the derivative 
operator 

d ^aodt + a^dz ^ ( ^0 9 " ) ' 

with d"^ = {n/c)dt ^ dz. can be thought of as basis vectors for a 

tangent space TAi of the two-dimensional manifold A4 formed by the (z, ct)- 
plane equipped with a space-time dependent refractive index function n{z^t) 
\W\ . The spatiotemporal evolution of the dielectric structure is only contained 
in the logarithmic term 9 Inn, which couples the forward- and the backward- 
propagating components of the electromagnetic field. All conclusions of this 
paper are based on an analysis of the main Eq.([2]). 

One can appreciate the remarkably high level of symmetry between space 
and time contained in Eq.([2]) with the following arguments. The propagation 
of the local field occurs only along the (local) direction specified by 
while the field F~ evolves only along d~ . Therefore and d~ determine, at 
each point of the (z, ct)-plane, the local light cone along which the forward and 
backward components of the electric field propagate. The geometrical meaning 
of our construction based on Eq.([2|) is depicted in Fig.[Tl where examples of the 
local light cones defined by the directions determined by and the propagating 
forward/backward fields are shown (see also caption for further explanation 
of the geometrical meaning of our construction). 

Contrary to a widespread approach [17J, we do not make use of any complex 
slowly varying variable, such as a field envelope, and therefore Eq.([2|) remains 
valid for all regimes of propagation, and all quantities appearing in Eq.([2]) are 
real. A slow variable version of Eq.([2]) will be examined in Section [H and a 
comparison between the exact approach of Eq.([2|) and the approximate one of 
Eq. (|45|) will be presented. 

Note that the coupling between and F~ is only induced via the logarith- 
mic term 9 Inn in Eq.([2|), so that only space-time variations of the refractive 
index are responsible for coupling forward and backward waves during propaga- 
tion. We can therefore identify two basic cases in which the analytical solution 
of Eq.Q is straightforward, namely the case where n changes abruptly from ni 
to n2 along a line in the (z, ct)-plane, and the case of a homogeneous and static 
medium, where Eq.([2|) simply becomes 9^ = 0. In the following two subsections 
these two cases are examined individually. 

2.2 Calculation of the Interface Matrix 

We assume that the refractive index changes in a step-like fashion along a non- 
stationary interface in the (z, ct)-plane, postponing the discussion of the validity 
of this sharpness assumption to the end of the next paragraph. The generic 
form of the refractive index is then given by n(z, t) = ni + (n2 — ni)0(cos(^)z — 
sin(6>)ct), where ni,2 are the refractive indices of the first and the second media, 
and 9 is the Heaviside function. 6 denotes the angle between the ct-axis and 
the interface in the (z, ct)-plane, which we will occasionally also parameterize 
via the slope parameter /3 = tan(6>). Note that the quantity /3c represents the 
velocity of the nonstationary interface. This velocity is arbitrary, and in par- 
ticular it can be greater than the local speed of light. This is not in confiict 
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Figure 1: The geometry of propagation of fields along the directions deter- 
mined by is shown, where j = {1,2} refers to the refractive index nj. The 
grey area represents the spatiotemporal region of refractive index change in the 
plane (z^ct). For this simple example, the external and the internal regions 
have uniform indices ni and n2 respectively. Bold dashed arrows indicate the 
propagating fields, generated by the scattering of F^ from an interaction point 
on the interface. Local light cones are also indicated. Assuming that n2 > ni, 
wider light cones are associated with index ni, whereas narrower light cones are 
associated with index n2. Only space-time variations of refractive index are able 
to generate coupling between forward and backward waves in Eq.([2]), so that 
in the case shown in the figure only the interface is able to scatter an incident 
electromagnetic wave. 
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with special relativity, because the nonstationary interface is not associated with 
a moving medium, and it constitutes one of the main differences between the 
present work and previous literature on nonstationary dielectrics. For different 
aspects of the physics of moving media the reader is referred to chapter IX of 
Landau [4|, section 1.5 of Jackson [3], and to more recent works by Yeh and 
collaborators [5j, Saca [6j, Huang [7] and by Leonhardt and Piwnicki [9j. 

Our task is now to find an appropriate transfer matrix connecting the am- 
plitudes of the fields on either side of the interface. We call this an interface 
matrix. Using the appropriate orthogonal coordinates 



p{z^t) = cos{0)z — sm{0)ct = 
q{z^ t) = sin(^)z + cos{6)ct = 



(3ct 



-13^ 
ct 



(5) 
(6) 



normal (p) and parallel (q) to the interface in the (z, ct)-plane, we can express 



the derivative operator in Eq.([2]) as d 
defined as 



Adp + Bdq , with the matrices A and B 



A = [az cos{0) — n<Jo sin(^)] 
B = [cr^ sin(^) + ncTQ cos(^)]. 

Taking into account that dqU = we can write Eq.([2|) in the form 
dp + A-^Bdq + A-^(Ja:A{dp \nn)cT^] ^ = 0. 



(7) 



(8) 



Writing EqE] in the basis [nE,B]^ 
ordinary differential equations: 



{ax -\- az)"^ yields the following set of 



dp {\nnE) 



- [-/3{l^n^)dq{\nnE) 



daB 



{uE) 
1 



-n (1 + /?^) dq {nE) + 2np (dpu) (nE)] 



-[-/3{l^n')dq{B) 



(9) 



(10) 



0. 



We can now integrate Eq. ([9]) from p = —e to p = +e across the interface 
at p = 0, where e is an arbitrarily small positive number. If B and {nE) 
are differentiable with respect to q (i.e. dqB < oo and dqlnnE < oo), and 
additionally ^ 1 across the interface, then the contributions to this integral 
from terms containing dq{B) and dq (nE) are proportional to e and will therefore 
vanish in the limit e ^ 0; Eq. ([9|) can then be integrated analytically. Under 
the same conditions Eq. (p!Q|) can also be integrated analytically and we obtain 
a linear relation between [niEi, Bi]^ and [n2£^2,^2]^- In the basis 
this relation can be expressed via the interface matrix D^^^] for the step-like 
transition that connects the fields F^ in medium 1 with those in medium 2: 



£)[12] 



F,^ 



(11) 
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Explicitly we find 




(12) 



This interface matrix possesses the classical property {D^^"^^) ^ = D^^^]. For 
future reference it is also useful to calculate 

detZ)l-l = ri^)!^. (13) 

Note that in the derivation of [12] we have used that n^/?^ 7^ 1 as we cross the 
interface. In particular this implies that our derivation is not valid in the regime 
where l/n2 < \P\ < l/ni. The physical reason for this will be discussed at the 
end of the next subsection, after we introduce the free propagation matrix. 

2.3 Free Propagation in Real Space 

From the knowledge of the general transfer matrix for a nonstationary interface 
[Eq. (p!2|) ]. it is possible to construct the transfer matrix for more complicated 
nonstationary objects, such as multilayer structures, which will be done in Sec- 
tion [3l However, we still need to construct one missing ingredient, the so-called 
propagation matrix^ which models the propagation of the fields along the p- 
direction far from the interfaces, where the refractive index is constant. Under 
this condition, Eq.([2|) simply becomes 9^ = 0, which using the definitions of 
Eq. [71 can be expressed as 

Adp^ ^Bdq^ = 0. (14) 

To solve Eq. ([T4|) . we first write the derivatives dp^q as first order finite differences, 
namely = [^(p+Ap, g)]/Ap and = [^{p, q^Aq)-'^{p, q)]/ Aq, 

and then substitute into the equation: 

^^ ^{p^Ap,q)-^{p,q) ^ ^{p,q^Aq)-^{p,q) 

Ap Aq ' K ) 

For a fixed Ap, we can eliminate the term proportional to ^(p, q) in Eq. ([T5|) with 
an appropriate choice of Ag, different for the two components of the ^-vector. 
From Eqslllone obtains 

^ I sin(^) zb ncos(^) ^ , 
cos(6') ^ nsm(6') 

where the refractive index n refers to the homogeneous layer under consider- 
ation, and Aq^ refers to the component F^. With the particular choice of 
Eq. ffTH]) . the free propagation of the optical fields for constant refractive index 
assumes a particularly simple and effective form, which we write for the two 
components: 

F^{p+Ap,q) = F^{p,q + Aq^). (17) 

This means that advancing the field along the p-direction by a small step 
Ap (keeping q constant), is equivalent to advancing the same field by an amount 
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equal to Aq^ along the g-direction (keeping p constant). The geometrical in- 
terpretation of Eq. (pT|) is that the points (p + Ap, q) and (p, q + Aq^) belong to 
the same light cone. 

The concept of the interface matrix [Eq. (p!2|) ] expressed in Section 12.21 is 
only valid for a sufficiently sharp interface, i.e. if the interface transition layer 
is sufficiently thin, such that the effect of the propagator (pT|) can be neglected 
while traversing the interface. For a given interface width Ap we can calculate 
via (p!6|) the length of the interval A^^ along the interface which is relevant for 
the propagator. If does not change significantly on this interval, i.e. if 



: ± n 
dq\nF^-—-—Ap 



« 1, (18) 



then the propagator is negligible and the use of the interface matrix of Eq. (p!2|) 
is justified. Note however that the left hand side of EqilHl becomes singular for 
np = ±1, and then [TSl will not be fulfilled for any choice of Ap. In particular, 
for n2 > ni, if l/n2 < /3 < 1/ni, then during the integration of [8] we will 
reach the point where the singularity in [18] occurs. This is precisely at the point 
where the matrix A~'^ also becomes singular. The condition on the sharpness 
of the interface which was used in the derivation of Eq. (p!2|) can therefore in 
principle not be fulfilled and therefore we shall not use Eq. (p!2|) in the regime of 
l/n2 < (3 < 1/ni. 

2.4 Plane wave expansion. Calculation of Propagation 
Matrix 

So far we have expressed the interface matrix of Eq. (p!2|) and the free propagator 
of Eq. (pT|) as operating on the fields vector ^(p, g), where p and q are the rotated 
variables of Eqs. ([5]) and (|6]). However, the interface matrix of Eq. (p!2|) connects 
the fields and F~ only along the p-direction, while not affecting the fields 
along the g-direction at all. It is therefore natural to consider a plane wave 
expansion along the g-axis, 

*(p,q) = ^*-(p)e--'^. (19) 

a) 

Then Eq. (pT]) transforms into 

^u^jj[12]^u^ (20) 

which means that the generalized frequency oj is conserved across the inter- 
face. Co has the dimension of a wavenumber, and physically it represents the 
momentum associated to the propagation along the (^-direction. 

Inside a homogeneous medium we can expand ^(p, g) in plane waves along 
both the p and q axes: 

^(p, q) = Y^ ^^rk^iCkp-^q). (21) 

oj,k 

Since each plane wave component has to fulfill the propagator equation d^^^^e^^^^ 
0, we obtain for the two components F^ the following dispersion relations: 

l± _ sin(l9)±ncos(l9) ^^ 

^ " cos(^)Tnsin(^)'^' ^^^^ 
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which can also be directly deduced from Eq.([T6]). This enables us to rewrite the 
free propagator (pT|) for the Fourier components as 

vl/-(p + Ap) = PiAp)^^ip) = j f -(p), (23) 

an expression that will be useful in Section [321 

As a final issue, it is worth writing explicitly the transformation law between 
the momenta reciprocal to the real-space variables (p, q) respectively, and 

(/c, cc;), reciprocal to the variables (z, ct). By using Eqs.([5l)-(|6l) and the expansion 
of Eq. ([2T]) we obtain: 

f ~k\_f cos(^) sin(^) \fk\ 

\Co J ~ \ -sin(^) cos(^) J [co/c J' ^^^^ 

From the simple rotation of Eq. ([24|) one can deduce some important physical 
consequences. First of all for 6 = {P = 0) it follows that k = k and uj = u/c, 
so that in the static case the p-direction coincides with the z-axis and the q- 
direction coincides with the ct-axis by definition. But because the g-direction 
is associated with the delocalization direction of an arbitrary plane wave in 
Eq. ([2T]) . it follows that during the scattering process of light by a static interface 
the quantity uj = uj/c is conserved. This corresponds to the energy conservation 
at the interaction point. In other words, all generated waves will have exactly 
the same frequency. The other limiting case is when ^ 7r/2 {/3 ^ oo). For this 
case the role of frequency and wavenumber is reversed, and from Eq. ([24|) we have 
k = uj/c and uj = —k. The delocalization direction is still along which now 
points along the z-axis. Therefore during the scattering it is the wavenumber 
(i.e. the momentum) that is conserved at the interaction point, and in general 
all generated waves will have different frequencies but the same momentum. For 
an intermediate situation (0 < < 7r/2), the conserved quantity Co corresponds 
to a combination of k and co/c according to the rotation of Eq. ([24|) . 



3 Applications 

In the present section we discuss three fundamental applications of the concepts 
expressed in Section 12. 2[ and in particular of the transfer interface matrix given 
by Eq. ([12]) and the free propagator given in real space by Eq. (fTTj) and in Fourier 
space by EQ. (f23|) . 

In subsection 13.11 we analyze the physics of scattering by nonstationary in- 
terfaces in different interesting physical configurations, shown in Fig. [2l In 
subsection 13.21 we calculate the bandstructure and the forbidden bandgaps of 
a nonstationary photonic crystal, finding a novel effect: the closure of the fre- 
quency bandgap for a certain value of ^ = 6>cri, and the consequent opening of 
the wavenumber bandgap for another = 6cr2- In subsection 13.31 we show that 
it is possible to think of a dynamical device, analogous to a spatiotemporal lens, 
that efficiently performs pulse broadening and compression, with simultaneous 
frequency /wavenumber manipulation of pulses. 
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0<|P|<l/n2 
(d) 













K n, 








l/n,<|P|<oo 



Figure 2: Panel showing space-time diagrams of scattering of forward and back- 
ward fields (j = {1,2}) from the four most important interface config- 
urations, (a) Space-like interface (/3 = 0, ^ = 0), which coincides with the 
conventional, static interface, (b) Interface configuration with /3 < and 
< \P\ < l/n2. (c) Time-like interface, P = ±oo {0 = ±7r/2). (d) Interface 
configuration with P < and 1/ni < < oo. 
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% 4 

Figure 3: (Color online) Panel showing the simulation of interaction of an incom- 
ing pulse with the four sharp interface configurations of Fig. O (a) Space-like 
interface: = 0. ih) [3 = —0.3. Angle = —0.2915 radians is indicated with 
dashed lines, (c) Time-like interface: (5 = oo. id) (5 = —2. Angle = —1.1071 
radians is indicated with dashed lines. Refractive indices are: ni = 1 and n2 = 3. 
The agreement with theoretical expressions is excellent for all cases. Input pulse 
parameter w in Eq. (|33|) is w = 1. Horizontal axis is dimensionless space coordi- 
nate ^ = z/A, while the vertical axis is dimensionless time coordinate r = ct/A, 
and A is the central pulse wavelength. 
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3.1 Nonstationary interfaces 

The panel of Fig. [2] shows the four most important cases of interface configu- 
rations which can be analyzed with the help of Eq.([T2]). Let us first consider 
Fig. [2](a) which shows the space-time diagram of a conventional static interface 
= 0, /3 = 0) between two media of refractive indices ni and n2 respectively in 
the plane (z, ct). Because all the variation of the refractive index is spatial, and 
the variable p defined in Eq. ([5|) can be identified with the spatial coordinate 
we call this a space-like interface. By taking the limit P ^ in Eq. (p!2|) we find 
the well-known result for the interface matrix: 

The transmission and refiection coefficients are then classically defined for F2 = 
as 

* - ^ = -^ = -^. (26) 



F2^ 


1 


2n2 




^11 


ni + n2 


Fr 


£,[12] 


ni - n2 


F,^ 


£,[12] 
-^11 


ni + n2 



(27) 

For < 1/3 1 < l/n2 as in Fig. [2](b) the configuration is similar to the one in 
Fig. [2](a). From the interface matrix Eq. (p!2|) we can now calculate the gen- 
eralized transmission and refiection coefficients according to the definitions of 
Eqs.([26D-([2lD as 

^ _ 2n2 1 - nip ^ _ ni - n2 1 - nip 
ni + 712 1 - n2P ' ni + 722 1 + niP ' 

Note that for < and inside the parameter region < \P\ < l/n2, the 
expressions given by Eq. ([28|) never show any singularity in the denominator. 
However, for > 0, we have a singularity in the transmission coefficient for 
P 1/^2- This limit corresponds to the disappearance from the equations of 
one of the waves participating in the scattering (^2^), which will travel exactly 
along the interface, so that t as given by Eq. ([28|) cannot be used, and for P > 
l/n2 the incident light never reaches the interface. 

For \p\ =00 Fig. [2](c) shows the configuration when the refractive index 
changes abruptly and simultaneously in time throughout all space (for some 
discussion of the experimental relevance of this concept see the conclusions. 
Section [5]). Because the variation of the refractive index occurs only along the 
ct-axis, we shall call this a time-like interface. From Eq. (p!2|) we obtain the 
interface matrix 

^[12] ^J_fni^n2, n2-ni\ 

2niV^2-ni, ni^n2 )' ^ ^ 

For the calculation of the transmission and refiection coefficients we have now 
to set = 0, and then instead of Eqs. (f26|) and ([27|) we obtain 

t ^ 4 = M'/' = ^^, (30) 
F+ " 2n2 ^ ^ 

, ^ ^=42i]^!^l^. (31) 
F+ 2n2 ^ ' 
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Similar formulas for this special case have been given by Morgenthaler in 1958 
[T2]. For a more general \/3\ > 1/ni when the interface is tilted with respect 
to the time-like case, the configuration is depicted in Fig. [2](d). Then Eqs. ([3Q|) 
and (|3T]) generalize to 

1 - pni ni + 712 _ 1- pnin2- ni 
~ l-(3n2 2ni ' ^ ~ 1 + (3n2 2ni * ^ ^ 

Even in this case, the expressions of t and r in Eq. ([32|) have singularities for 
(5 = l/n2 and (5 = — l/n2 respectively, which however are outside the region of 
validity of these formulas, that is \(3\ > l/ni. 

In Section 12.41 we have seen that the frequency uj arising in the plane wave 
expansion along q in Eq. ([2Q|) is conserved across the interface. As we have 
already discussed previously, in the case of /3 = [Fig.[2ja)] this is, with a) = cj/c, 
the well-know frequency (i.e. energy) conservation at a spatial interface. In the 
more interesting case P = oo [Fig. EJc)], however, we see that Co = — /c, which 
means that wavenumber (i.e. momentum) is conserved at a time-like interface. 

The scattering of finite pulses from the interfaces for all four cases shown in 
Fig. [2] have been simulated by solving numerically Eq.([2]), and results are shown 
in the panel of Fig. [3l Electric field for the input pulse is taken of the form 

Eini = E{z, t = 0)=Eo sech cos (^^^ , (33) 

where Eq is the incident amplitude, and w is di dimensionless parameter which 
measures the spatial width in units of A, the pulse central wavelength. Without 
loss of generality, the incident amplitude is in all cases normalized to the value 
^0 = 1. 

We have used a second order finite difference algorithm to model the deriva- 
tives in Eq.([2|) for the two components, and a fourth order Runge-Kutta algo- 
rithm to advance the field in time. Once the spatial distribution of the electric 
field over the z-axis at the initial instant of time Eini = E{z,t = 0) is given, 
the initial magnetic field Bini = B{z^t = 0) can be found by imposing that 
there is no backward wave propagating for t = 0, i.e. F~ = 0, obtaining 
Bini{z) = n{z^t = 0)Eini{z)^ where n{z^t = 0) is the initial refractive index 
distribution. We use Neumann-type boundary conditions, which means that we 
impose that the derivatives of the optical fields vanish at the boundaries of the 
spatial window, while the evolution of the fields is calculated by advancing in 
time. In order to ensure complete stability of the numerical algorithm it is suffi- 
cient to impose that all fields are spatially localized and that the temporal step 
of the grid At and the spatial step of the grid Az satisfy the following condition: 
At <C Az'^. For our simulations we always use a spatial sampling of the optical 
fields with a number of points equal to = 2^^, which is more than sufficient 
in order to obtain extremely accurate results for long propagation times, and to 
avoid the dangers of artifacts due to numerical dispersion during the temporal 
evolution. 

The panel of Fig. [3] shows propagation of the envelope of the electric field 
(with the fast oscillations decoupled for sake of clarity) during the scattering 
processes analogous to Figs. [2fa,b,c,d) respectively. Here and in the following, 
contour plots are given in the dimensionless variables r = ct/X and ^ = z/A, with 
A = 27rc/a;o, and coq is the central circular frequency of the initial pulse. Note 



13 




Figure 4: (Color online) (a) Scattering diagram in the case < 0, l/n2 < < 
1/ni, and (b) corresponding numerical simulation. Parameters are: /3 = —0.4 
{0 = —0.3805 radians, indicated with dashed lines), w = 1. Refractive indices 
are: ni = 1 and n2 = 3. 

that Eq.([2]) is solved directly in dimensionless units, therefore we do not need to 
specify the scaling wavelength A, and results shown in the figures can be adapted 
to the particular spatial scale that one chooses to consider. The extension of 
the spatial grid used is (^window = 300 in dimensionless units, which is relatively 
large in order to ensure total spatial localization of all waves participating in 
the scattering, although only a small portion of this window is shown in the 
figures. The spatiotemporal steps are = 9.15 x 10~^ and Ar = 7.67 x 10~^. 
The pulse is initially propagating at the speed of light in vacuum (ni = 1). The 
sharp interface is modeled through a Heaviside function n{z^t) = ni + (n2 — 
ni)Q{z — /3ct)^ with n2 = 3. The values of /3 for Fig. [2](a,b,c,d) are respectively 
P = 0^ P = —0.3, = oo, and P = —2. In all cases, t and r as given by the 
theoretical analysis given above are predicted with a high degree of accuracy. 

As discussed at the end of Section [231 the derivation of the interface matrix 
(p!2|) does not cover the case of l/n2 < \P\ < 1/ni. The configuration for this 
case is depicted in the upper panel of Fig. [H We see that this case resembles the 
space-like configuration as in Fig. [2ja) in the region n = ni and the time-like 
configuration as in Fig.[2jc) for n = n2. Note that now only one of the four fields 
{F^) reaches the interaction point from the past, while the other three fields 
propagate from the interface in positive time direction. It is therefore not clear 
how in this case the transmission or refiection coefficient can be consistently 
defined, and neither of the given definitions of Eqs. ([26|) - ([27|) nor Eqs. ([3Q|) - ([3T]) 
are satisfactory in this case. The theoretical origin of this problem is that in 
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the region l/n2 < \/3\ < 1/ni the transfer matrix concept is not apphcable, 
due to the fact that for these values of /3 the concept of sharp interface cannot 
be defined, at least using a dispersionless susceptibility, which does not make 
distinctions between long and short wavelengths. In particular, the condition of 
Eq. (p!8|) does not necessarily hold in the above region, as recognized by Ostrovsky 
[llj. Nevertheless we can still solve Eq.(|2]) numerically as shown in Fig. [H^b), 
and we leave the investigation of this complicated issue for a future publication. 



3.2 Nonstationary Photonic Crystals and Bandgap Cal- 
culation. Spatiotemporal Bragg reflector 

We now use the knowledge built in Section 12.41 to calculate the bandstructure 
and the forbidden bandgaps of a nonstationary multilayer structure for different 
values of and in particular we shall enunciate the Bragg condition for a 
generalized Bragg r effector. 

Let us assume in the following analysis that the periodicity of the multilayer 
structure is along the p-direction [see Fig. [5ja) and caption for a description 
of the geometry of the problem under consideration] . As in Section 12.41 we now 
deal with the propagation of plane waves of the form ^f{p^q) = ^^{p)e~^^^. 
Similarly to the extensively studied Kronig-Penney model [18 , we now look for 
a matrix T such that 

^'^(p + Ap) =T^'^(p) (34) 

which propagates ^^{p) from the center of a layer with refractive index ni 
across a layer with refractive index n2 into the center of the next layer with 
refractive index ni as indicated in Fig.[5l[a). The length of this unit cell is given 
by Ap = pi + p2 , where pj are the widths of the layers of indices Uj along the 
p-direction. From Eqs. (f2Q|) and ([23|) we obtain 



T = Pi(pi/2)Z)[i21p2(p2)Z)[2ilPi(pi/2) (35) 

where Pj{p) is the propagation matrix according to Eq. ([23|) in the medium 
with refractive index rij. Since Dt^^l^t^^] = 1 we observe that |detT| = 1. 
Considering the two eigenvalues Af and \2 of T we can then distinguish between 
the following two cases: (i) |Af | = |A2^| = 1 and (ii) |A|^| > 1 and |A2^| < 1. 

Free propagation will only occur in case (i), which by elementary algebra 
can be shown to be equivalent to 

i|Tr(r)|<l. (36) 

In this case it is convenient to write the eigenvalues of T in the form 

Ai = e^'k^ve^^K A2 = e^~^^ve-^^K^ (37) 



where k is calculated via VdetT = e'^^p and K fulfills the relation 

Physically k simply gives rise to a common phase factor in both eigenvalues 
of T, while K can be interpreted as a generalization of the Bloch wavenumber 
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in the Kronig-Penney model for nonstationary interfaces. In particular the 
explicit calculation of K via a straightforward analytic expansion of the trace 
and determinant of T yields the familiar- looking expression [18j: 



cos(i^Ap) 



= cos(piA;f^) cos(p2^2^) 



(39) 



sin(pifcf ) sin(p2fc#), 



where 




cos- 



■?{e)-sm^{0)n] 



rijuj 



(1 + l3^)njCb 
l-/32n2. 



(40) 



For a given Co we can now use Eqs. (|39l) - (|4Q]) to calculate the appropriate 
if it exists. This gives rise to a band structure relation as shown in Fig. [6] (see 
also caption). Note that [3 only enters in the band structure calculation via k^^ 
which however diverges for j3 l/n^ [cf. Eq. (|4Q|) ]. 

The forbidden bandgaps are located in those regions of Co for which condition 
(|36|) is not fulfilled, i.e. the right hand side of Eq. (|4Q|) is larger than unity and K 
would have a non- vanishing imaginary part, making the Bloch waves evanescent. 
In particular if the quantities pj are chosen such that 



then one can form a Bragg refiector with total refiectivity inside the forbidden 
bandgaps in Co. Condition (|4T]) is analogous to the quarter wavelength condition 
typical of a Bragg stack. However, and this is the crucial point, it should be 
noted that the meaning of Co is different for different values of (5. For instance, in 
the traditional case of space-like Bragg stack, i.e. in the limit ^ 0, condition 
(jUJ) becomes the well-known uodjUj/c = 7r/2, because in this limit rijCo^ 
Co ^ uojc and dj (the width of layer j), where uo is the initial plane wave 

frequency. However, in the time-like Bragg stack, i.e. in the limit (3 ^ oo [see 
Fig. [51(b)], Eq. (|4T]) becomes ktj/rij = 7r/2, because in this limit kj^ Co/rij, 
Co ^ —k^ and pj ctj (the duration of layer j) where k is the initial plane wave 
wave number, common to all waves. 

Fig. [6] shows the bandstructure (solid lines) and the forbidden bandgaps 
(grey areas) in Co for a generalized Bragg refiector that satisfies the generalized 
quarter wavelength condition of Eq. (|4T]) . with 0<^<7r/2, as calculated by 
solving numerically the transcendental photonic band equation Eq. ([39|) . For 
this special case, all bandgaps have exactly the same extension and the spacing 
between them is regular [2 . All momenta in Fig. [6] are normalized to dimension- 
less units by multiplying by Ap, and K is the Bloch wavenumber. Orthogonal 
basis (a), k) is rotated with respect to the basis (a;/c, k) of an angle ^, indicated 
in the figure, see also Eq. ([24|) . If ^ ^ 0, a) ^ cj/c, the two axes coincide and 
the bandgaps in Co would correspond exactly to the bandgaps in uo/c. This is of 
course the space-like case of a static, conventional Bragg refiector. For the other 
extreme case, ^ 7r/2, Co ^ — /c, and the bandgaps in Co would correspond to 
bandgaps in k. This novel configuration corresponds to a time-like Bragg re- 
fiector, of the type shown in Fig. [H^b). In the intermediate cases < ^ < 7r/2, 
depicted in Fig. [6l the projection of one of the bandgaps onto the co axis is 
shown, and its extension is marked with a dashed double-arrow. It is clear that 
for increasing values of the extension of the projected cj-bandgap is shrinking. 



Pjkf = 7r/2, 



(41) 
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Figure 5: (a) Schematic figure showing the quantities of interest used in the 
text to define a general multilayer structure (with indices n^, j = {1,2}) for 
an arbitrary angle 6 (i.e. arbitrary parameter P). p indicates the propagation 
direction for the Bloch waves, while q is the delocalization direction for plane 
waves, dj are the widths of the layers along the 2;-direction, while ctj are the 
widths along the ct-direction. pj indicate the same quantities but along the p- 
direction. Ap = pi -\- p2 is the total length of the unit cell, also indicated in the 
figure, (b) Time-like photonic crystal. The refractive index varies only along 
the time axis, which therefore coincides with p. The coordinate z coincides with 
therefore plane waves are spatially delocalized. For the structure to be a 
time-like Bragg reflector, which is totally reflective for a certain wavenumber 
bandgap, the condition ctj = rijX/A must hold, where A is the wavelength of 
the incident plane wave. 



until for a certain critical angle Ocri its size vanishes, indicating the closure of the 
frequency gap. A further increase of will lead to a second critical angle ^cr2, 
in correspondence of which the bandgaps projected onto the /c-axis gradually 
start to open, until they reach their full extension at ^ = 7r/2. 

With the above discussions we have therefore identified the novel general 
concept of a spatiotemporal Bragg reflector^ and hence of a spatiotemporal pho- 
tonic crystal with a generalized quarter wavelength condition given by Eq. (|4T]) . 

3.3 Spatiotemporal Lenses, Pulse Compression and Broad- 
ening 

Another interesting and potentially useful effect exhibited by a simple class 
of spatiotemporal dielectric structures is what we have called spatiotemporal 
lensing. Let us suppose that we have a refractive index distribution in the 
(2;, ct)-plane of the form: 

n(z, t) = ni + (n2 — ni) exp 

In Eq. (|42|) . zq and to represent the coordinates of the center of the rectangular 
region, and Wz^t are respectively its spatial and temporal extension, see Fig. [71 
The two integer numbers m^,^ model the sharpness of the supergaussian transi- 
tions between ni and n2 along the spatial and temporal directions respectively. 
Eq. (|42|) models a rectangular region in the (z, ct)-space where the transition 
between an 'external' refractive index ni and an 'internal' index 77,2 occurs, see 



2mt 



Wt/2 J 



2mz 



(42) 
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Figure 6: Bandstructure (solid lines) and bandgaps (grey areas) of a spatiotem- 
poral Bragg reflector that satisfies the generalized quarter wavelength condition 
of Eq. (|4T]) . as calculated by using the transcendental photonic band equation 
Eq. (|39|) . According to Eq. ([24|) . the orthogonal axes lj and k are rotated with 
respect to the orthogonal axes u and k by an angle 6>, as indicated in the fig- 
ure. All quantities shown are dimensionless, and normalized in units of Ap. For 
^ = (space-like photonic crystal), the two orthogonal bases would coincide, i.e. 
Co = co/c and k = k^ and this corresponds to the case of a conventional, static 
photonic crystal. For a general 6>, the forbidden bandgaps in lj as calculated by 
Eq. (|39|) are shown by the grey areas. The projection of one of the bandgaps 
onto the axis co is also shown, with its new smaller extension indicated by a 
dashed double-arrow. In the limit ^ 7t/2 (time-like photonic crystal), uj tends 
to — /c, while k tends to cj/c, so that the forbidden gaps in Co will correspond 
to bandgaps in k. Note that there is a critical angle = Ocri at which the 
projected cc;-gap will close, and a second critical angle 6cr2 > ^cri at which the 
/c-bandgaps start to open, until they reach their full extension at 6> = 7r/2. 
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Figure 7: Schematic figure showing the mechanism of pulse compression [(a)] 
and pulse broadening [(b)] due to spatiotemporal lensing. The refractive index 
of the rectangular region shown is given by Eq. (|42]) . The refractive indices 
outside the rectangle and inside it are uniform, in the (z, ct)-plane, and equal 
to ni and n2 respectively. In this example we assume that n2 > ni. During 
the scattering process every portion of the pulse will travel according to its 
local light cone, which is narrower in the grey region (light travels at speed 
c/n2 < c/ni), and wider outside it. Backward waves produced in the scattering 
are not indicated in the figure for sake of clarity, but they are always present 
for sharp interfaces. 

Fig. [71 Here and in the following we shall refer to this region as a spatiotemporal 
lens, for reasons that will be clear shortly. 

Fig. [71 illustrates the essential features of the scattering of a generic pulse 
by the rectangular spatiotemporal lens in the (z, ct)-plane. Essentially two con- 
figuration are of interest. In the first configuration [Fig. [71(a)], after some free 
propagation in the medium with index ni, the pulse enters the grey region with 
index n2 on the left side of the rectangular region, with the entire body of the 
pulse well inside the lateral extension of the rectangle (wt). The pulse then 
propagates inside the structure, at a greater angle in the {z, ct)-plane, each por- 
tion of the pulse following the local light cone (defined in Fig. [T]) at any instant 
of propagation. When the pulse exits the grey region, however, it is clear from 
Fig. [T^a) that its pulse duration has been reduced, and pulse compression has 
been achieved. It can be shown by elementary geometrical considerations that 
the compression ratio between the output and the input pulses in this simple 
case is just p = ni/n2. The scaling in the (z, ct)-plane is accompanied by an 
increase of the central frequency of the pulse by with a consequent increase 
of the spectral bandwidth. Thus we obtain an effective method for manipulating 
the frequency and wavenumber of a pulse. 

The second scenario is when the pulse is delayed in such a way that it enters 
the spatiotemporal structure on the bottom side of the rectangle, the extension 
of which is and well inside it, see Fig. [71(b). This case shows exactly the 
opposite dynamics of the previous case, and pulse broadening in the (z, ct)- 
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Figure 8: (Color online) Simulation showing pulse propagation and scattering 
by a spatiotemporal lens in the (z, ct)-plane which separates two regions with 
uniform refractive indices ni and n2, as schematically described by Fig. [71 (a) A 
broad pulse [w = 3 in Eq.(j33l)] hits the rectangle on the leftmost side, and pulse 
compression is achieved, (b) A short pulse [w = 1 in Eg. ([33]) ] hits the structure 
on the bottom side, and pulse broadening is achieved. For sake of clarity, only 
the interface of the rectangle is shown. Parameters of the simulations are: 
ni = 1, n2 = 3, zq = 0, cto = 20A, Wt — 27A, = 15A, rrit = rriz = 20, 
and A is an arbitrary scale which we take equal to the central pulse wavelength. 
Horizontal axis is dimensionless space coordinate ^ = 2:/A, while the vertical axis 
is dimensionless time coordinate r = ct/A, and A is the central pulse wavelength. 

plane occurs by a factor of p = 712/^1. As a consequence the central frequency 
is reduced. This is of course not surprising due to the symmetry of Maxwell's 
equations under the time reversal operation. 

In Fig. [71 the backreflections at the interfaces, when the pulse enters the 
structure and when it leaves it, are not shown for simplicity, but it is clear from 
the discussions of the previous sections that for sharp interfaces of the structure 
they are always present. 

The scenarios described above are confirmed by direct numerical simulations, 
see the panel of Fig. [H Fig. [81(a) shows the evolution of a spatially and tempo- 
rally broad input pulse [w = 3 in Eq. (|33l) ]. After hitting the leftmost side of the 
rectangular spatiotemporal region, the pulse emerges from the upper side con- 
siderably compressed, both in spatial and temporal domain. Only the interface 
of the 712 region is indicated for sake of clarity, in order to clearly distinguish the 
internal wave propagation. Fig.[8l^b) shows the compression dynamics, in which 
an input pulse with w = 1 hits the structure on the bottom side. The pulse 
emerging from the right-hand side acquires a significant broadening. Several 
waves propagating backwards are noticeable in Fig. [8l^a,b), and these are due 
to the sharpness of the interface region, which is modeled by Eq. (|42l) with the 
specific parameters given in the caption. 

So far we have considered the extreme cases where both interfaces are space- 
or time-like. With the formalism developed in the previous section we can also 
consider the general case where a pulse enters a medium with refractive index 
712 via an interface characterized by an angle Oi and leaves that medium via 
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Figure 9: (Color online) Comparison between pulse evolution described by (a) 
the full Eq.([2]) and (b) the approximate Eq. (|45|) . which was derived under as- 
sumption of SVEA in time. The two simulations are carried out under absolutely 
identical conditions. Parameters are: P = —0.2, = —0.2027, = 1, ni = 1 
and 712 = 3. In (b) one can notice an unphysical artificial broadening along the 
z-direction due to the uncompensated term proportional to the second deriva- 
tive in space 9| in Eq. (|45|) . Horizontal axis is dimensionless space coordinate 
^ = z/A, while the vertical axis is dimensionless time coordinate r = ct/A, and 
A is the central pulse wavelength. 



another interface characterized by an angle ^2- A straightforward geometrical 
analysis gives the rescaling factor in the (2:, ct)-plane as 



cos c/2 — ni cos Ui 
cos O2 — n2 cos Oi 



sm c/2 — n2 sm Ui 
sin ^2 — ^1 sin^i 



(43) 



which as expected reduces to p = ni/n2 ( p = n2/ni ) for the previously 
considered cases Oi = and O2 = 7r/2 {O2 = and Oi = 7r/2). 



Derivation of a Slow- Variable equation. Com- 
parison with Eq.([2]) 



From Eq.([T]) it is possible to derive an equation for the electric field envelope 
function A{z, t) only, which is of second order in the variables z and t, by 
posing E{z^t) = [Aex.p{ikoz — iuoot) + c.c.]/2, where ko = ujq/c and ujq is an 
arbitrary frequency, usually and most conveniently chosen to be close to the 
central frequency of the input pulse. Using the dimensionless variables r = ct/A, 
^ = z/A, where A = 27rc/ct;o, and substituting the expression of the electric field 
E into Eq.([T]), we obtain: 



ndrU 



drA-'^^idlA) 



■ id^A + —d'iA+ 



2indrn {druf nd^n 



TT 



27r2 



27r2 



47r 



A = 0. 



(44) 
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Although Eq. (|44|) is written for the slow variable A, it is exacts and it is able 
to handle arbitrary pulses, and forward and backward propagation at the same 
time, due to the fact that it is completely equivalent to the full wave equation 
which can be derived from combining the two Maxwell equations of Eq. ([T]) . 

It is reasonable to conjecture that the use of the slowly varying envelope 
approximation (SVEA) in time in Eq. (|44|) is valuable in order to reduce the 
computational complexity of the problem. However, as we will show shortly, 
the use of SVEA in time but not in space is dangerous in that it leads to wrong 
results for pulse propagation, and this is especially important for the dielectric 
systems we are investigating in this paper, which show strong variations in the 
refractive index both in space and time. The use of SVEA in space in Eq. (|44|) 
is not permitted in particular when simulating sharp interfaces and photonic 
crystals, and in general in those systems in which the backward propagating 
waves are of the same order of magnitude of the forward propagating ones. 
Therefore all second derivatives 9| must be maintained in Eq. (|44|) . 

Assuming that the central frequency of the pulse ujq is much greater than 
the variation of the envelope function A in time, we can perform a multiple scale 
reduction of Eq.(j44j), obtaining the following first order equation in time: 



A - 

47r ^ 



1 



2indrn 



A = 0. 



(45) 



Eq. (|45|) is a generalization to time-dependent refractive index of the propagation 
equation used in a series of works [17J, where temporal variations of n were not 
taken into account. Eq. (|45|) shows that the purely real 'Schroedinger potential' 

— 1 used in [17J is modified by a purely imaginary term proportional to the 
time derivative of the refractive index, namely idr'n'^ /tt. 

Fig. [9] shows the comparison between the scattering of a pulse by a sharp 
interface, when calculated by means of the full Eqs.([2|) [Fig. ^a)] and by means 
of the approximate equation Eq. (|45|) [see Fig. [H^b)], under absolutely identical 
conditions. Although the two equations give very similar results for what con- 
cerns the transmission and reflection coeflicients, the most striking difference 
between the pulse evolution in the two cases is the presence of a fictitious and 
unphysical broadening along the 2;-direction in Fig. [9jb). This is due to the fact 
that the use of SVEA in time but not SVEA in space in Eq. (|45|) leads to an un- 
compensated 'diffr active' term proportional to which is responsible for the 
broadening observed in Fig. [91(b) for all waves participating in the scattering. 

The above considerations provide a justification of the use of the full Maxwell 
equations throughout the present paper, instead of equations using slow vari- 
ables. 



5 Discussion and Conclusions 

Before concluding, we would like to spend a few words on the experimental 
relevance of the topics discussed in the present paper. 

Many theoretical and experimental results exist in the literature, that dis- 
cuss the interaction of electromagnetic pulses with a plasma front [T9l [201 [2H [22 } 
[23l [24j . This interaction has been typically used in practice to shift the central 
frequency of an incident pulse by means of the effect discussed in Section [331 
see Refs. [21 , or to convert static electric fields into electromagnetic pulses by 
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means of a capacitor array [23] . These effects may turn out to be extremely 
useful in the creation of efficient tunable laser sources For instance, in typ- 
ical experiments [2^, when a 'probe' pulse is propagating along the longitudinal 
coordinate of a semiconductor waveguide, one can change the refractive index 
of the material along all its spatial extension by means of irradiating the plane 
of the crystal with another 'source' pulse that creates a electron-hole plasma 
throughout the crystal. If the source pulse has a certain angle of incidence a 
with respect to the normal to the plane of the waveguide, then the ionization 
front of the plasma will move with velocity v = c/ sin(a) p^^J. Note that for 
normal incidence (a = 0), this front moves instantaneously, and it corresponds 
to the time- like interface {0 7r/2, where the general relation between and 
a is -\- a = 7r/2) discussed in Section ISTTl which follows from a simultaneous 
change of refractive index throughout all the spatial extension of the crystal. 
This front does not carry any information, and it is not associated to any moving 
object, therefore its propagation is not limited by the speed of light in vacuum. 
In all these cases, none of the actual particles in the excited medium is moving 
faster than c; instead, the front profile only has an apparent motion, which can 
be characterized by an effective velocity given by the above formula, but does 
not correspond by any means to any real superluminal motion of particles or 
information, which is of course forbidden by special relativity. In this sense, 
n{z,t) can be a totally arbitrary function with no relation whatsoever to mov- 
ing bodies. In another kind of experiments [24], short pulses of the order of 
nanoseconds are created by irradiating a metal surface by X-rays at a certain 
angle of incidence a, which creates a photoelectric current. Again, this current 
undergoes an apparent motion at a velocity v = c/ sin(a) and, by means of su- 
perluminal self-phasing of the excited elementary dipoles (predicted by Carron 
and Longmire in 1976 [24J), emits a short electromagnetic pulse at the same 
angle a. It is also interesting to point out that the processes described briefly 
above can also have considerable implications for the physics of astrophysical 
plasmas [23]. 

A remark which is important for the present work is that it is assumed that 
the physical experiments well approximate a condition of quasi-one-dimensional 
propagation of the incident probe pulse, therefore implying a well-localized con- 
finement of the transverse profile by means of a strong waveguiding process 
which does not allow diffraction, so that the equations we have derived in this 
paper, having the transverse (x, y) and the longitudinal {z) degrees of freedom 
decoupled, can be safely utilized. 

In conclusion, in this paper the physics of light propagation scattered by 
nonstationary interfaces has been investigated analytically and numerically in a 
unified approach. Interface transfer matrix and propagation matrix have been 
found for the general case of arbitrary velocities of the interface. These ingre- 
dients have been used to construct more complicated spatiotemporal dielectric 
structures. In particular, we have ffist investigated a spatiotemporal Bragg 
reflector, where we have given the generalized condition for the existence of 
forbidden bandgaps, and found that the frequency bandgaps close for a critical 
value of ^ = ^cri, with the successive opening of the wavenumber bandgap at 
a second critical angle = 0cr2' Secondly we have explored light scattering by 
a spatiotemporal lens device, and the spectral manipulation of pulses has been 
demonstrated both analytically and numerically, the most important aspect of 
which is a potentially strong pulse compression/broadening. 
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Future works include the natural extension of the theory given here to a 
dispersive susceptibility. 
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